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1. Introduction Ls ction ASL: t ret Fa pe tik Ie a Aa 1% 
., Lanchester~type* differential~equation combat models are.an important 


tool for, analyzing many important problems of military operations research. 
In such a combat model, a so-called attrition-rate coefficient represents” 
the fire effectiveness of a particular weapon-system type against a particular 
target type, i.e. the weapon-system type's effective firepower against such 
a target. Time-dependent attrition-rate coefficients are used to model! 
temporal variations in firepower on the battlefield. Thus, we see that 
time~dependent attrition-rate coefficients are importent (and, in fact, 
essential [4-6)) for the quantitative analysis of hypothetical combat. 
Militarily realistic computer-based Lanchester-type models of quite 
complex military systems have been developed for almost the entire spéctrum 
of combat operations, from combat between battalion-sized units to theater- 
level operations. Nevertheless, a simple combat model may yield a clearer 
understanding of significant interrelationships that are difficult to per- 
ceive in a more complex model, and such insights can subsequently provide 
valuable guidance for more detailed computerized investigations. In this 


report w2 consider such a simplified variable-coefficient Lanchester-type 


' model of combat between two homogeneous forces. 


For thie variable-coefficient Lanchester-type model of combat 
between two homogeneous forces, different functional forms for the attrition- 
rate coefficients lead to different mathematical functions being involved 
in representing and computing the force-level trajectories. In a previous 


_ paper {5] we have discussed the plausibility of the hypothesis that except 


for the special case of a constant ratio of attrition-rate coefficients, 


* so-called after pioneering work of F. W. Lanchester [3]. 


the solutions to such differential equations cannot be represented in term 
of "elementary" functions of analysis. Thus, new transcendental functions — 
arise in the study of combat modelled with time-dependent attrition-rate 


coefficients. In particular, we have previously introduced’ [5-6]' so-called 


Lanchester-Clifford-Schlafli (LCS) functions for analyzing combat modelled’ 
ithe power attrition-rate coefficients with "no offset" (see Section 3 below). 
In the Appendix to this report is contained a reduced set of tables 

for the LCS functions: it contains tables of five-decimal=place values of 

the hyperbolic-like LCS functions F(x), H, _, &), and T(x) (see Section 4 
below) for 11 fractional values of a (see Section 6 below). A companion 
report [8] contains the most extensive set of tables currently available. 

The main body of this report provides the theoretical and modelling background 
for the use of these tables. In particular, we examine a model of a constant- 
speed attack ona static defensive position and show how associated range- 
dependent kill rates give rise to time-dependent attrition-rate coefficients 
with "no offset." Numerical computations are presented to illustrate the 

use of the LCS functions for analyzing such "aimed-fire" combat. As a con- 
sequence of the availability of these tables, one can now study this variable~ 
coefficient. combat model almost as easily and thoroughly as Lanchester's 


Classic constant-coefficient model. . 


2. Variable-Coefficient Lanchester-Type Equations of Modern Warfare. 


We consider combat between two homogeneous forces modelled by the 


following variable-coefficient Lanchester-type [3] (see [4,5]}) equations 


of modern warfare 


NAS eR - ae Dee at i i Ni a 


dx 


7 * ~a(t)y “with x(0) = Xo » 
: (2.1) 
Lm (t)x with (0) = y5 


where t = 0 denotes the time at which the battle begins, x(t) and y(t) 
denote the numbers of X and Y at time t, and a(t) and b(t) denote 
time-dependent Lanchester attrition-rate coefficients, which represent the 
sitecitecunes of each side's fire. These coefficients depend on variables 
such as force separation, tactical posture of targets, rate of target acqui- 
sition, firing rate, etc. (see [4-7] for ticther detatis). Variable attrition- 
rate coefficients are used to model temporal variations in firepower on the 
battlefield. In any analysis of combat, moreover, we should use the above 
equations (2.1) only for x and y > 0 and, for example, set dx/dt = 0 
when x = 0, since negative force levels have no physical meaning. 

Mathematically, we assume that the attrition-rate coefficients a(t) 
and b(t) are defined, positive, and continuous for ty <t<+eo with 
ty < 0. We also assume that a(t) and b(t) € L(t),T) for any finite 


T2 tp. We further take a(t) and b(t) to be given in the form 
a(t) = k g(t) , and b(t) = kh(t) , (2.2) 


where k, and k, are positive constants chosen so that a(t)/b(t) = k,/k, 


when g(t) = h(t). We introduce the combat-intensity parameter A, and 


the relative-fire-effectiveness parameter Ap defined by 


a a Se da tae ers i 


pe? 


, of time is given by a similar expression, with C(t) and S(t) being 


4, © vic k, : and  * ki /k, (2.3) 


From our assumptions about a(t) and b(t), it follows that, for example, -4 


T 
a(t) f L(ty,T) implies i a(t)dt = +0. 
0 


The X force level as a function of time may be represented 
as [5,6] 
x(t) = xp{C) (OC (t) - S,(O)S,(t)} - yoVAQ {C,(O)S,(t) - S.(O)CL(t)}, (2.4) 


where the hyperbolic-like general Lanchester functions (GLF) Cy (t) and Sy (t) 


are linearly-independent solutions to the X force-level equation 


ax 
aa at 


ate: x - a(t)b(t)x = 0, (2.5) 
with initial conditions 

C, (to) cate Se S, (to) = 0 
(2.6) 


{1/a(to)} de/dt(t)) = 0, {1/a(t,)} dS /dt (ty) = Vr, ; 


Here ty denotes the largest finite time at which a(t) or b(t) ceases 


to be defined, positive, or continuous. The Y force level as a function 


analogously defined for the corresponding Y force-level equation. 


It is sometimes convenient to introduce the new independent variable 


t defined by 


» 
t= f{ Va(s)b(s) ds . (2.7) 
*9 
It is readily seen that the transformation 1 = 1t(t) is well defined and 
invertible. Let us denote 1(0) as Tp: We observe that ty <0 implies 
9 2 9- If we denote the “average intensity of combat" as 
Ya(t)b(t) , then 


that T 


Fob at SO t ; 
Ya(t)b(t) t = {(1/t) { Va(s)b(s) ds}t =1- Typ (2.8) 
0 


The substitution (2.7) transforms (2.5) into 


2 
ax _ Afa Ee zh 
e- {Ss tn R(t)} SP - x= 0, (2.9) 
dt 
with initial conditions 
x(t) = Xp >» and {1/VR(xp)} dx/dt(t)) a a's 


where R(t) = a(t)/b(t). 


3. Combat Modelled with Power Attrition-Rate Coefficients. 
me EES with Fower Attrition-Rate Coefficients 
The above equations (2.1) basically apply to “aimed-fire" combat 


when target-acquisition times do not depend on the numbers of targets available 


(see [5,6] for further details). A large class of tactical situations of 


interest can be modelled with the following general power attrition-rate 


coefficients [5-7] 


a(t) = k(t + c)", and b(t) =k (e+ C+ a. G.1) 


where A and C > 0. ‘We will call A_ the offset parameter, since it allows 
us to model (with uy and v > 0) battles between opposing weapon systems 
with different maximum effective ranges (see [5,6]). We will call C the 
starting parameter, since it allows us to model (again, with wu and v > 0) 
battles that begin within the maximum effective ranges of the two opposing 
systems. We observe that for the general power attrition-rate coefficients 


(3.1) we have t, = -C, and uw and v must be > -1 in order that a(t) 


0 
and b(t) € L(t),T). 
The above nomenclature is motivated and possible applications of our 


work are indicated by considering S. Bonder's model of the constant~speed 


attack on a static defensive position (see [4-7] for further details) 


ax. 2 ae 

at a(r)y, and ae B(r)x , (3.2) 
where xr denotes the range between opposing forces, and a(r) and 8(r) 
denote range-dependent attrition-rate coefficients. Range is related to time 
by 


r(t) =R, -vt, (3.3) 


0 


where Ro denotes the opening range of battle and v > 0 denotes the 


_constant attack speed. For example, let us consider the constant-speed attack 


of a homogeneous Y force against the static defensive position of a homo- 


geneous X force. Figure 1 diagrammatically portrays this situation. 
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The basic idea is that force separation, i.e. range between the opposing 


forces, changes over time, and the fire effectiveness of, for example, 


a single Y firer, denoted as a(r), depends on this force separation. 
In many cases of tactical interest, we may model the fire effective- 


ness of, for example, the Y weapon system (as a function of range) with 


q r lg ‘ 
7 | ay 5 ee a for O<r<R, 
| a(r) = (3.4) 


0 for R. <r, 
a 


where R denotes the maximum effective range of the Y weapon system and 
u > 0. Here yp is used to model the range dependency of Y's attrition-rate 
coefficient (see Figure 2). We model B(r) similarly, with corresponding 
quantities Re and v being analogous to R and u_ above, 


If we use (3.3) to eliminate range r from (3.4), we obtain 


TE = -a(t)y , 


; (3.5) 


dy 2 
at b(t)x , 


where the time-dependent attrition-rate coefficients a(t) and b(t) are 


given by (3.1). It follows that the offset and starting parameters are given 


‘by 
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R, -R R -R 
As ( Pent 1, and C= (4.-%o), (3.6) 


and that : ee | 


u ; sid ss v . 
k, = 45 (z) » and ky = By (x) p (3.7) | 
a B 


We observe that A and C > 0 4f and only if R, > a2 R,. By consider~ 


8 0 
ing (3.6) and Figure 3, the reader should have no trouble in understanding 
our terminology for A and C, 

When the offset parameter is equal to zero (i.c. A = 0), then the 


iy 


coefficients (3.1) reduce.to 


a(t) = kj(t+0)", and Bt) = W(t +0)” (3.8) 


We will refer to (3.8) as power attrition-rate coefficients with "no offset." 


As we have seen above in Bonder's constant-speed attack model, these 
coefficients model, for example, combat between weapon systems with the same 
maximum effective range so that there is no "offset" in the "reaching out" 
of the weapon systems against each other in combat (again, see Figure .3). 
For these coefficients (3.8), the transformed X force-level srecten (2.9) 


becomes 


2 
dx 2q-1, dx ica g 7 
me 2 + ¢ . ) ae 0, (3.9) 


with initial conditions 
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) ce ‘ c) + v+2)/2 


, a 
3 - 1, + 3)e (ut v+2)/2 ; (3.12) 


‘ Let us observe ‘that 0 <q<il \when u and v>-1. Furthermore, q > 1/2 
- if and only if “ aR/at <0, i.e. R(t) is a strictly decreasing function of 
» time. fae 


Consider the function Fi (x) di fined by the power series 
“ 8 t ‘ Me baa 
\ si 


ae SS oo : 2k ™ 
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For a ¥ 0, °-1, m2reh the radius of convergence for F,(x) is in 
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| by the ratty test for clavicgiacage power ~~ [2]. — F, (2) te 
Q 
| an entire finction of the souphen Meciiide z= = a+ ty with an essential 
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singularity at the point at infinity. Now consider the function H, (x) 


ia ane HR ea 


defined by the infinite series 


= xa i 2(kta) i | 
: H (x) = T(a) J ea : (4.2) 
F a k=O k! Tik+a +1) | 
: | 
7 Observing that | 
; Hox) = (1/o) (x/2)7" FL Gx) (4.3) 
; a atl . é 
7 we see that for a >0O the infinite series (4.2) is uniformly convergent 
: on compact subsets of the complex plane. From (4.3) one can readily deduce 
3 . | the recursive relation 
F(a) = F(x) + ol Woy (4.4) 
a atl [ata +1) { “a+2 : : 
We will call the functions F(x) ‘and H (x) Lanchester-Clifford-Schlafli 
(LCS) functions (see Note 10 on pp. 66-67 of [5]). Other properties are 
readily deduced and are given in Table I. 
The function F(x) satisfies the linear second-order ordinary 
differential equation ae Sas Eh oy pt 
i 2 
| dF dF 
, —% 4 @el 8p . 
F a oto ae = 7, 2 8; (4.5) 


a with initial conditions 
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Table I. Properties of the LCS Functions P(x) and H, (x). 


1-2 
1. dF /dx = (x/2) : “Hy 


2a-1 
2. dH /dx = (x/2)" “F, (x) 


cA reenter 


3. Fy(x)F,_g(x)-H,GB,_ Gx) <2 Vx 


where a is not an integer (including zero) 
Fe 4. Fg Gr =1 
‘cbe Hy (x0) = 0 for a>Q0O 
6. dF /dx(x=0) = 0 


1-2a 


Te {G0/2)" "an fax) oy = 2 


8. Fy sa = cosh x 


9. By 72 (*) = sinh x 


ar, : 
F,,(0) = 1, and = (0-0, 


while BO) satisfies 


di 
- - Ch gt-a-o, —" 


with initial conditions 


te a, 


ax =l. 


x=0 


BA(0)=0, and 16 


Thus, {F., By.) is a fundamental system of solutions to 


ere a=}) = -F=0, (4.7) 


with Wronskian W(F_, H, |) = (x/2)*""*_ rt follows that the CLF for the 


X and Y force-level equations for combat modelled with the attrition-rate 
coefficients (3.8) are given by 


. ee 2q-1 
I . 
C(t) = F(r(t)), 5, (t) o_o } H(r(e)) ’ (4.8) 


' mt £ 


' | ee 1-2q Ju 
Cy(t) = F (reed), S(t) 7 B (v(t) » (4.9) 


. 


where p= lq. If we define 


T, (x) = Hy ng 2) / FY (x) , (4.10) 


then 
s ot) -( Ay me H (r(t)) 
0° ay Soe ees | FG@) Ty) ada 
or 
ry ae “ 
T, (t) = rare sy E Ty ¢r(e)) ‘ (4.12) 


where T,(t) denotes a hyperbolic-like GLF, which corresponds to the 


hyperbolic tangent. Observing that for wu, v > -l, lim T(t) = +, 


TH +o 


we see that TG) is a strictly increasing function of x on the interval 


{0, +) and 
o<t,() < for O<x<to, (4.13) 
with 
lim) T(x) = 2G) 4.14 
din T0025} (4.18) 


since by the results of Taylor and Comstock [7] the parity-condition parameter 


Q* = Qk(y, v, C #0) is given by © | ; 


d 2q-1 
eee 
oo a Oe o - ( ) | (4.15) 


We recall that Taylor and Comstock [7] have introduced the so- 
called parity~condition parameter Q* as the value (or range of such values) 


for the initial condition Q to the initial-value problem 


Bai A it Rk, tt tats 


: dE, 1 e 
5 (4.16) 
dE inn es . 
— =~ dp b(t)E, with Ey (t)) =Q, 


such that EY (t3Q*) and EY (t3Q*) > 0 for all t >t). In other words, 

Q* is the value of Q in (4.16) above such that neither EY nor Ey 

ever become zero. In this case, both EY (tsQ*) and EY (t3Q*) are positive, 
strictly decreasing functions, similar to decreasing exponentials. Thus, 

we may call Q* "the Y equivalent of an X force of unit strength," 

since the forces are “at parity," with neither force being annihilated in 
finite time, Taylor and Comstock have shown that for either a(t) ¢ L(0,+ ~) 


or b(t) £ L(0, + ~), then Q* is unique and given by 


Sy (t) iL 
lim F (4.17) 
te+o e C(t) * Qt 
The significance of the parity-condition parameter Q* is that it allows 


us to predict force annihilation as the following theorem shows. 


THEOREM 1 (Taylor and Comstock [7]): Assume that either 
a(t) £ L(O, +”) or b(t) & L(O, +”). Then the X force 


will be annihilated in finite time if and only if 


Cy (0) - Q*s, (0) 
. (4.18) 


*o AAR 
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Use of LCS Functions for Analyzing Combat. 
The Lanchester-Clifford-Schlafli (LCS) functions PG) and H, (x) i 


i ie ii 


are useful for analyzing "“aimed-fire" combat (see Section 3 above) modelled 


i 
| 5. 
| 
: 
with the power attrition-rate coefficients with "no offset" (3.8), which 


we rewrite here as 


prensa ines 7 nret 


a(t) = k(t +0)", aaa b(t) = k(t +0)”. (5.1) 


In other words, the LCS functions arise in solving the differential combat i 


model (2.1) with attrition-rate coefficients (5.1). Im order that both 


a(t) and b(t) € L(tp,T), we must have wy and v >-l. Military situations 
modelled by these equations have been discussed in Section 3 above, e.g. 
combat between two weapon systems with the same maximum effective range. 
For such combat, the LCS functions may be used to 

(1) compute force-level declines, 

(2) predict force annihilation, 
and (3) predict the time of force annihilation. 

Let us now see how the LCS functions may be used to obtain the 
above information about force-level declines and force-annihilation pre- 


diction. According to (2.4), (4.8), and (4.9) above, the xX force level 
is given by : 


x(t) = Xo fF, (t) F (x(t) ~ H(t) A (r(t))) 


I 2q-1 
- Yo rp <spet (F(t) H(t (t)) HC) Fo (r(t))}, (5.2) 
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where q is given by (3.10), p = 1-q, and 7+(t) is given by (3.11), which 


we rewrite as 


2d 
i I (u4v+2)/2 
T(t) i .) (t + C) F (5.3) 


The time to annihilate the X force* is determined by x(t%) = 0, and 


it follows that 


i Ay yes 
oF“ + ¥% Malievea) Mp 


a q-P 
I 
Xgl (tT? + Yo vx (5) FAO) 


x 
Fi (r(t.)) = 


R\ptvut2 


where from (4.10) 
Ty Cr(t)) = H(t (t))/F C(t) > (5.5) 


and we recall that p+q#=1. It follows that the time to annihilate X, 


ef, is given by 


"Te we multiply the first equation of (2.1) by y, the second by x, add, 
and integrate, we obtain 


t 
x(t) y(t) = xpyq - { {als) y"(s) + b(s) x(s)}ds , 


* which shows that x(t) and y(t) can have at most one finite zero. Hence, 
ee x(tX) = 0, then we know that y(t) > 0 for all t > 0. 


QT 


r q-P 
I 
a x,F (t,) + y gy a H (t.) 


oY . 
I 
Xp (to) + ¥9 Mp (_ Paes Ui FAS? 


Taylor and Comstock [7] have shown that r&) is strictly increasing and 


satisfies (see also (4.12) above) 
0< tT) < I(p)/T(q) , (5.7) 


where p = J-q. It follows that in order for X to be annihilated in finite 
time, the right-hand side of (5.4) must be less than [(p)/I(q). Let us 


observe that for to = -C = 0, (5.4) simplifies to 


P-q 
x Xo Ay ) 
fe )) + ae 5 (5.8) 
se Yo Ye 


Thus, we have proved the following theorem concerning force- 


annihilation prediction. 


THEOREM 2: Consider combat between two homogeneous forces 
modelled by (2.1) with power attrition-rate coefficients 
(5.1). Assume that wp and v>-1 and that the above 


equations hold for all time. Then the X force will be 


annihilated in finite time if and only if 


My ye t 
PD) 1 Xo FoR Nieves! Hho) | 


TpLanhr EI MO WOOT ORG a 


q-P 
{ ( » 
< T'(p) | XoHg{o? + Yo Wap ieee S F (to) { , (5.9) 


where q = (v+1)/(u+v+2) and p= 1-q. For ty = 0 


(i.e. C= 0 so that Tp = 0), X will be annihilated in finite 


time if and only if 


Xo r(p) ( ry a 
Yo Fda) Wy \ eee ; (5.10) 


6. Tabulation of LCS Functions. 
This report contains a reduced set of tables of the Lanchester-Clifford- 
Schlafli functions. The Appendix contains tables of five-decimal-place values 
of the hyperbolic-like LCS functions F(x), Hy _,{*), and T, (x) for various 
values of the argument x, namely x = 0.00 (0.01) 2.00 (0.1) 10.6, and 
a = 1/2, 1/3, 2/3, 1/4, 3/4, 1/5, 2/5, 3/5, 4/5, 3/7, and 4/7. These values 
of the index a correspond to u, v = 0, 1, 2, and 3 in (3.8) and allow one 
to analyze, for example, a basic spectrum of range capabilities for weapon 
systems in the constant-speed-attack model of Section 3. These tables have 
Bk, been calculated by the recursive means given in Section 8 of [5]. A more 
| extensive tabulation (namely, for a = 1/2, 1/3, 2/3, 1/4, 3/4, 1/5, 2/5, 


3/5, 4/5, 2/7, 3/7, 4/7, 5/7, 4/9, 5/9, 3/11, 5/11, 6/11, 8/11, 5/13, 8/13, 


1 
? 


5/17, 12/17, 5/21, and 16/21 corresponding to u,v = 0, 1/4, 1/2, 1, 1 2, 3) 
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is to be found in a companion report [8]. This companion report contains 
the most extensive set of tables of the Lanchester-Clifford-Schlafli 


functions currently available, 


A representative tabulation of the hyperbolic-like LCS functions 
F(x), H,_,@&) and T(x) is given in, for example, Tables 8A and 8B of 
the Appendix for a= 3/5. The values of the argument x are the same as those 
used for the tabulation of the hyperbolic functions by Abramowitz and Stegun 
[1]. We observe from Table 8B and (4.13) that the limiting value of 
T, @&) as x*>+o (here a = 3/5) is quickly reached, with three-decimal- 
place accuracy already attained for x= 4.5. Moreover, the use of these 
tables (specifically, Tables 8A and 8B of the Appendix) for combat analysis 


is illustrated in the next section. 


7. Numerical Examples 


In this section we examine a couple of numerical examples to show 
some of the insights that may be gained into the dynamics of combat between 
two homogeneous forces from our results (see also [6]). These examples 
illustrate the use of the LCS functions F(x), H,_,@), and T(x) for 
analyzing "“aimed-fire" combat modelled with the power attrition-rate coefficients 
with "no offset" (5.1). As in [4-7], we consider S. Bonder's model (3.2) 


for the constant-speed attack against a static defensive position. We will 


focus on the use of the LCS functions for predicting force annihilation, 


since the computing of force-level trajectories with Lanchester functions 


is adequately handled elsewhere (see [4~5]), 
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Let us accordingly consider the constant-speed attack of a homogeneous 
Y force against the static defensive position of a homogeneous X force 
(see Section 3 above for further modelling details, especially Figure 1). 
For our numerical computations, we assume that the fire effectiveness of the 


Y weapon system varies linearly with range, i.e. 


r 
9 (1 - x) for Sse sh, 
a(r) = : (7.1) 
0 fcr R, it > 


and that the fire effectiveness of the X weapon system varies quadratically 


with range, i.e. 


2 
ty (1 -) for O<r<R,, 


B(r) = (7.2) 
0 for R, € Big 
with R. = Ry» i.e. both weapon systems have the same maximum effective range. 
In other words, » = 1 in (3.4) and v=2 for 8(r). We consider a battle 
modelled by the input data given in Table II. In terms of time as the 
independent variable, the attrition-rate coefficients (7.1) and (7.2) become 


via (3.3) 


a(t) = k(t + ¢) and b(t) = k(t + er. (7.3) 


23 


RN ih NE Als ol Saeki a pA Be hte Sl Re PE a aati 


Table II. Input Data for Numerical Examples , ‘ 
| u® 1, v=2 
ay = 0.06 X casualties/minute/Y firer 
' 
| By = 0.6 Y casualties/minute/X firer 
| 
R =R, = 2000 meters 
j 
v = 5 miles/hour 
| 
| 
| 24 


where Ry = Ra, 


R. - Ry . Gov ~ 2 
ate ae ° ea ose ° and k s(x) ° (7.4) 


From the input data given in Table II, we compute the parameter values shown 

in Table III, since the transformed X force-level equation is given by (3.9) 
with q= (v+1)/(u+v+ 2), p = l-q, » = 1, and v= 2. Thus, the xX force 
level may be computed with F(t) and Hy _,() with a=q = 3/5. Force- 
annihilation prediction involves the limiting value of T,() = Hy _ (D/P (7) 
as t++, From Table 8B of the Appendix and Table III, we note the predicted 
agreement between [(1-0)/f(a) and the limiting value of T, x) as x*+t+° 
[recall (4.13)] for a= q = 3/5. We now sdaxtser two cases (1) Ry = 2000 
meters, and (II) Ry = 1250 meters. 


When R, = 2000 meters (see Figure 3 of [4]), we have C=O and 


0 
To = 0. The maximum time that the battle can last is t... = R/V = 14,91 minutes, 
since at this time the attackers reach their final objective, i.e. the 
defender's position (again, see Figure 1). We now consider the qualitative 
behavior of the » = 1, v= 2 force-level trajectory shown in Figure 3 of 
[4]. Theorem 2 tells us that the X force can be annihilated in finite time 


1f and only if 


Xo Tp) Ay q~-P 
Yo PCa) fi (aetiey) , (7.3) 


where q = 3/5 and P = l-q. Using the numerical values in Table III, we 


compute from (7.3) that the X force can be annihilated in finite time if 


‘and only if 


“or 
4g , 
i 
ej 
a! 
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Table III. Parameter Values for Numerical Examples 


= 4.0233 x 107? x ‘casualties/(minute)"fy firer 


o5% ‘ ? a ih eye <3 } < be geeny By Fe os Ege iry Bes see 
Va mays Ss & Vip Siete >. BESS > wes A ests eS pd eS y 34 > Bars i Be 


“ky = 2.6979 « 107 ¥ casualties/(minute)”/X firer 


yey eats z -¥ vat } me tf Ce : Ltat Tes SRS BE 


% 2S 4 $6 SGuphee ertb FWatsl aes ssviowice wot seipe ste  Shise hie 
p= 2/5, q = 3/5 
ge OS eS ae ee Pa Sore eee? eee bie tT EEE Heows et Ape tas 
- P(p)/T(q) = 1.48951 
rag: "Ge : 7 ae i a WS eS iS ie aes j ft LtF 4 Rae | 
A = 0 4 Sa E i " Pi t2 is eras 3 


ee arc 


x 
5, < 0-420. (7.4) 


¥ 


When the X force can be annihilated, its annihilation time is given by 


(5.8), which we rewrite here as 


x a P-q 
T (r(t®)) = —2 (ts) : (7.5) 
q a Yo a utvt2 
where 
2d 
I (utvt2) /2 ‘ 
t(t) = eee Gog ‘i (7.6) 


Thus, for the numerical values given in Table III, the time of annihilation 


of the X force is given by 


x a) 
Ty (r(tQ)) = 3.544 7” (7.7) 
"y 


with q= 3/5. We will now illustrate further computations for Xp» * 10. . and 
Yo * 30. From (7.4) we see that the X force can be annihilated in finite 
time (but we must verify that * £ tax) * In this case (7.7) becomes 


¥ 
Tyre) = 1.18122 . fe ke (7,8) 


We must now determine 1(t*) such that 1(t™) = 7 (1.18122) by using inter- 
polation methods and Tables 8A and 8B. From Table 8A, we find 
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Tho = 1.18172 for = 1,01 


TT) = 1.17630 for = 1.00 
so that using linear interpolation, we obtain 
e ; 
T(t) = 1.009 , (7.9) 


whence use of (7.6) yields 


t® = 14.24 minutes, (7.10) 
which is less than ne = 14.91 minutes so that he defending X force is 
indeed annihilated before the attacking Y force reaches its final objective. 

: | Since r(t) = Ry - vt, we find that force separation at the instant of hs 


annihilation of the xX force is 
xX . 
* 89.8 meters . - (7.11) 
Further results may be computed in a similar fashion and are given in 


Table IV. 
When R, = 1250 meters (see Figure 3 of [5]), we have C = 5.5923 


minutes, Tt)" 0.0975, and fas = o/¥ = 9,32 — In this case 
Theorem 2 tells us that the X force can be annihilated in finite time 


P| ‘4 and only if 
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Table IV. Annihilation of the X Force as a Function 


of the Initial Force Ratio for Ro = 2000 meters 


(x9/Y¥ 9) (minutes) rs (meters) 


0.333 0° 16.24 89.8 
0.250 11.61 443.2 
0.200 10.19 633.2 


PAR te Pore erermereeartereene verter. 


rey cy} 


utvet2 T(q) Pocftenor , 
‘ p‘to) T'(q) H(t») } 


with q = 3/5 and p= 1-q. Using linear interpolation, we obtain from 


Tables 7A and 8A of the Appendix that for the numerical values of Table III 
F(t) = 1.006 , H(t) = 0.044 , 
F(t) = 1,004 , H(t) = 0.223 , 


so that (7.12) says that the X force can be annihilated if and only if 


When the X force can be annihilated, its annihilation time is given by 


(5.4), which we rewrite here as 


*o ( Ay 
utv+t2 
Yo VR 


P-4q % 
Fy ft) ae 2 (. a +2 Hy (9) 
Yo"'R 


P-q 
F(t)? + ss 


x 
T (r(es)) = 


whence for the data of Table III 


x. 3+565uy + 0.223 
T,(t(t,)) = 0.156, + 1.008” 


aici accata RAN A BAS 8k ark F AE a eA IRS Li eaitaip lin W am AOA WII aS 


where Up * Xo/¥o: Let us also record here that (3.11) yields 


2/ (utvt2) 
t= (fep+2) na (7.17) 


I 


We will again illustrate further computations for X> * 10 and Yo * 30. 
From (7.14) we see that the X force can be annihilated in finite time (but 
again we must investigate whether or not ss < tax)" In this case (7.16) 


becomes 
7, (x(t2)) = 1.33651 , (7.18) 
whence Table 8A of the Appendix and linear interpolation yield \ 
r(t2) = 1.397, (7.19) 
so that by (7.17) 
t* = 20.63 minutes . (7.20) 
Since t, | = R,/v = 9.32 minutes < at we see that the attacking Y force 


overruns the defender's position before annihilation of the X force occurs. 


Thus, the battle ends with X_ = x(t-) > 0 and Ye > 0 at a iets 


9.32 minutes. Corzvesponding to t. = 9.32 minutes is t= 1.1318, and 
then Table 8A of the Appendix yields 


Fa(t, = 1.1318) = 1.589 , H (1.1318) = 1.973 , 
whence via (2.4), (4.8), (4.9), and (7.13) we obtain 


x, = x(t ,) = x(r = 0) = 1.35. (7.22) 
Some further numerical results are given in Table V, Again, these parametric 
results should be contrasted with the single’ 1 = 1, v #2 force-level 


trajectory shown in Figure 3 of [5]. 


8, Final Remarks ‘ 
In the previous section above, we have seen how the LCS functions 
allow one to conveniently obtain much valuable information about the model 
(2.1) with power attrition-rate coefficients (3.8) without having to explicitly 
compute the entire force-level trajectories. Previously we were limited to 
computing only force-level trajectories (see [4-5]). With the availability 
of these tabulations of LCS functions (see the Appendix of this report and 
[8]), we can now tell who is going to be annihilated and when this event will 
happen without having to compute the trajectories. Not only did we answer 


questions about the qualitative behavior of the model (e.g. force annihilation) 


for specific values of, for example, initial force levels but also for a 


range of values of the initial force ratio (i.e. parametric analysis, of 


model behavior). 


Table V. Annihilation of the X Force as a Function 


of the Initial Force Ratio for Ro = 1250 meters 


(x9/Y 9) “4 (minutes) 


0.333 10.63 
0.250 7.56 


0.200 6.17 


+t 


on = 9,32 minutes and X_ = x(r=0) = 1.35. 


The results of this report may be used for other parametric 


analyses, e.g. parametric dependence of battle outcome on attrition-rate 
coefficients. Thus, ‘hk: desecaiaie of this report allow one to develop 
important insights into the dynamics of combat between two homogeneous 
forces with temporal variations in fire effectiveness. With the avail-~ 
ability of tabulations of the LCS functions, one can now analyze such 
combat modelled by the power attrition-rate coefficients (3.8) with some- 
what the same facility as he can for the constant-coefficient case and thus 
aid in parametric analyses. For further discussions of the significance 

of such results for military operations research, the reader is directed 


te [6) and [7]. 
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APPENDIX; Tabulation of the LCS Functions F(x), H,_, > and T(x) for 


a= 1/2, 1/3, 2/3, 1/4, 3/4, 1/5, 2/5, 3/5, 4/5, 3/7, and 4/7. 
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